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In order to increase the luminosity of electron-positron colliders it is desirable to 
find a means to reduce the phase-space of the beams. The transverse cooling of 
positrons imposed by the quantum excitation-free radiation emission in a single 
crystal is considered as a potential route to achieving ultra-cold beams. An anal- 
ysis of the problem is presented, including an evaluation of the contribution from 
multiple scattering during the passage. The analysis shows that an omittance re- 
duction may be achieved in special cases, but in general the emittance will increase 
as a result of the multiple scattering. 



1. Introduction 

In a series of theoretical papers, the quantum excitation-free radiation emis- 
sion in a single crystal has been discussed 1,2 . It is shown that the transverse 
action - and thereby the emittance - decreases exponentially towards the 
minimum value h/2, corresponding to an emittance of half of the Comp- 
ton wavelength. This applies as long as the radiation is in the undulator 
regime where the angle of emission is larger than the pitch angle. On the 
other hand, experiments 3 have been performed which in agreement with 
theoretical expectations based on a completely different analysis show that 
positrons in contrast to electrons generally suffer heating instead of cool- 
ing. The present work is motivated by two things: The desire to estimate 
the potential of single crystals as 'devices' for the production of ultra-cold 
beams and a wish to find a consensus between two apparently different the- 
oretical approaches and experimental results. In the following the outline 
of 2 is followed, with the inclusion of multiple scattering, and it is shown 
that for the experiment one should indeed expect heating as concluded from 
both theories. Finally, it is shown that in special cases it is expected that 
transverse cooling can be achieved. 



1 
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2. Channeling, multiple scattering and dechanneling 

The large fields present near the nuclei in solid materials may in the case 
of single crystals add coherently such that a penetrating particle experi- 
ences a continuous field along its direction of motion - the so-called contin- 
uum approximation 4 . If further the particle is incident with a sufficiently 
small angle to a particular crystallographic direction, inside the so-called 
Lindhard angle, the negatively/positively charged particle is constrained to 
move near/far from the nuclei and the electron clouds surrounding these. 
This is the channeling phenomenon 4 which has found widespread applica- 
tions in physics. For a general introduction to channeling and applications, 
sec eg. 5 . The critical angle for planar channeling is given by 

V> P = ^AZ 1 Z 2 e 2 Nd p CaT ¥ /pv (1) 

where Z\e is the charge of the penetrating particle, Z^e that of the lattice 
nuclei, p the momentum, v the velocity, Nd p is the planar density of atoms, 
N being the atomic density and d p the planar spacing, C ~ V3 is Lindhard's 
constant and otf is the Thomas-Fermi screening distance. 
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Figure 1. The figure shows as a function of reduced transverse position, x/d p , a) the 
transverse potential in the Moliere and harmonic approximations and b) the electron 
density calculated from Poisson's equation. 
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In the continuum approximation, the resulting transverse potential leads 
to a continuously focusing environment in which photon emission may take 
place without recoil to the emitting particle, the recoil being absorbed by 
the lattice. This is the so-called 'semi-classical channeling radiation reac- 
tion' ' . In the previous papers on this phenomenon two main assumptions 
are made: The particle is moving in a harmonic potential and the energy of 
the photons emitted is small compared to the energy of the particle. Disre- 
garding for the moment the potentially important case of axially channeled 
positrons along a crystal axis, this leaves only planar channeled positrons 
since channeled electrons are in a strongly anharmonic potential. Axially 
channeled positrons may be confined to a region between certain strings if 
their transverse energy is very low, so-called proper channeled positrons. In 
this case the transverse potential can be well approximated by a harmonic 
potential, see e.g. 6 . As seen from figure la, the harmonic approximation 
is clearly well suited for positive particles. 
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Figure 2. The squared scattering angle per unit distance, A9 2 /Az, as a function of 
reduced transverse position, x/d p , for a p = 25 GeV/c beam in silicon. 



The transfer of channeled particles to states above the barrier, so-called 
random, is referred to as dechanneling. The length, Lb, over which a planar 
channeled beam of protons has been reduced to the fraction 1/e of the initial 
intensity by transfer to the random beam is given for 7 3> 1 by 7 , 8 : 

T _ 256 pv aTF^p 

D ~ 9^ 2 ln(2 7 mc 2 //)-l Z ie 2 U 

where I is the ionisation potential. Eq. (2) has been shown to be in good 
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agreement with measured values of Ld at room temperature 8 . Due to the 
dependence of electron density on transverse position, see figure la, the 
dechanneling process which is a result of the multiple scattering depends 
itself on the transverse position and therefore on the transverse energy. In 
figure 2 is shown the squared scattering angle per unit distance, A9 2 /Az, 
for a p — 25 GeV/c beam in silicon. For low transverse energies, the 
multiple scattering is dominated by the interaction with electrons and it is 
observed that A6» 2 /A z is almost a factor of two lower for the more accurate 
Moliere approximation than for the harmonic potential (for the Molierc 
approximation, see e.g. 9 ). 



3. Quantum excitation-free radiation emission including 
multiple scattering 

Following 2 the time evolution of the longitudinal Lorentz factor, j z = 
^/m 2 c 4 + p 2 c 2 /mc 2 = E z /mc 2 , in a continuous focusing environment is 
given by the differential equation 

= - TcGl r ut) (3) 

i.e. it couples to the transverse action, J x = E x /u z , which evolves as 

^ = -T c J x -\T cGl lJ 2 {t)J 2 At) (4) 

leading to two coupled differential equations, both in the absence of multiple 
scattering. Here r c = 2r e K/3mc is given by the focusing parameter, K, 
related to the transverse potential height, £/ , and planar spacing, d p , as 
J7o = K(d p /2) 2 /2 and G = Kjm?c 4 is a convenient constant expressing 
the focusing strength. 

The solution is given as 

Mt) - J x0 (l + ^7'o^o(l - exp(-r c t)))- 3 / 5 cxp(-r c i) (5) 



and 



5 



l z {t) - 7,o(l + W z0 6 2 {l - exp(-r c t)))- 4 / 5 (6) 



where 9 p o = \J 2GJ x o / 7^ q 2 is the initial pitch angle, J x o the transverse 
action and 7 z o the longitudinal Lorentz factor upon entry. 

To include multiple scattering we use the analysis for the dechanneling 
length, Ld, of positive particles with 7 » 1. The dechanneling process 
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arises due to a steady increase of the transverse energy imposed by mul- 
tiple scattering, i.e. the transverse energy increases as dE x /dz — Uq/L-d. 
Therefore the transverse action increases as dJ x /dt — Uoc/Lr>u) z , but since 
the effective dechanneling length depends on the transverse energy approx- 
imately as L\y = LdqUq/ J x u) z , the transverse action changes according to 

dJ x /dt = c/L m J x (7) 

Here Ldo denotes the dechanneling length for states where E x ~ Uq. This 
is found by dividing the squared scattering angle for x = by the average 
over x of the squared scattering angles shown in figure and multiplying the 
dechanneling length from eq. (2) by this ratio, i.e. 

/•dp/2 

Ldo = L D d p A6 2 (d p /2)/2 / A6 2 (x)dx (8) 

Jo 

Thus, by combining eqs. (4) and (7), the result is 

^ = - T C )J X - \v cGl ll\t)Jt{t) (9) 

whereas j z (t) remains unaffected. By a change of variables T' c — T c (l — 
c/T c Ldo and G = G(l — c/r c LDo) _1 the same type of coupled differential 
equations as eqs. (3) and (4) are obtained with solutions given by eqs. (5) 
and (6) with T' c and G instead of T c and G and with a new pitch angle, 
6' p0 = (1 — c/r c ZvDo)~ 1//2 $pO- To get the true pitch angle as a function of 
time the inverse transformation is applied, i.e. 6 P = (1 — c/T cLdo) 1 ^ 2 ^ 

where 6' p — \J 2G J' x /l^ 2 is the modified pitch angle. As expected and 
seen from the results below, a good measure of the depth at which the 
cooling effect starts to appear is given by ct c = c/T c . 



4. Results 

In figure 3 is shown the pitch angle calculated as a function of normalized 
penetration time, t/r c , for 25 GeV positrons in a (110) diamond, with 
and without inclusion of the multiple scattering. It is seen that while the 
cooling starts around t/r c = 0.1 without multiple scattering, it is postponed 
to values above t/r c ~ 5 when this additional effect is taken into account. 
Since there is always an incoherent contribution to the radiation emission 
which typically takes place as in an amorphous medium, it is important 
to note the scale of cr c compared to the amorphous radiation length, Xq, 
which for diamond is 122 mm. It is thus not possible to utilize planar 
channeling in diamond for cooling of a 25 GeV beam of positrons for an 
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angle of incidence near ip p . However, for smaller angles the cooling starts 
already at ~ 0.1t/r c , i.e. a 10 mm thick diamond would suffice to initiate 
the cooling. 
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Figure 3. The pitch angle in diamond (110) as a function of normalized penetration 
time, t/r c . The full drawn curve is calculated including multiple scattering while the 
dotted line excludes this contribution. The graphs a), b) and c) are for incidence angles 
ipp, tpp/2 and ij) p /&, respectively. 

In figure 4 is shown the different behaviours for several angles of in- 
cidence, calculated for tungsten (110). Strong cooling is found to appear 
early for small values of the angle of incidence, but for tungsten ct c ~ 8X 
which means that even for small angles there will be a strong influence from 
incoherent scattering. 

In figure 5 is shown results for silicon and germanium showing that for 
light materials, where the lattice is not very compact as in diamond, the 
influence of multiple scattering is much stronger than for heavier ones. 

To get an impression of the variation of the cooling effect as a function 
of energy in figure 6 is shown the results for 4 different energies. Clearly, 
as the energy increases, the influence of multiple scattering diminishes and 
cooling starts earlier. 
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Figure 4. The pitch angles in tungsten (110) as a function of normalized penetration 
time, t/r c , for 4 different angles of incidence as indicated. The positron energy is set to 
25 GeV and all curves are calculated including multiple scattering. 




Figure 5. The pitch angle in a) silicon and b) germanium (110) as a function of nor- 
malized penetration time, t/r c . In each case, the initial pitch angle equals the planar 
critical angle, ip p , and the full curves are calculated including multiple scattering while 
the dotted curves do not include this effect. 



As a crude estimate, the axial potential can be approximated by a har- 
monic potential with barrier height equal to that found by the Molicrc 
potential. The two kinds of motions in the planar and axial channeling 
cases are very different and it is not expected to get anything but an indi- 
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Figure 6. The pitch angle in tungsten (110) as a function of normalized penetration 
time, t/r c , for 4 different energies as indicated. In each case, the initial pitch angle equals 
the planar critical angle, tp p , and all curves are calculated including multiple scattering. 



cation of the magnitude of the effect from this analysis. However, in figure 
7 is shown the results obtained by this procedure for two energies and we 
note that in this case the characteristic cooling length is much shorter than 
Xq. This indicates that the much stronger axial fields may indeed provide 
substantial cooling. 
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Figure 7. The estimated axial pitch angle in diamond as a function of normalized 
penetration time, t/r c , for 2 different energies as indicated. In each case, the initial 
pitch angle equals the planar critical angle, tp p , and both curves are calculated including 
multiple scattering. 



Such an experiment has been performed by NA43 at CERN, see 3 , and 
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compared to a theoretical analysis by Kononets, see 10 . It was found, in 
agreement with this theory, that positrons suffered heating in contrast to 
electrons. In figure 8 is shown results under the same crude estimation 
method for the axial case as mentioned above. The actual value of the 
crystal thickness in the experiment is given by the vertical dashed line and 
it is seen - as observed in the experiment - that for all angles of incidence 
positrons suffer heating. 
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Figure 8. The estimated axial pitch angle in silicon as a function of normalized pene- 
tration time, t/r c , for 4 different angles of incidence. All curves are calculated including 
multiple scattering. 



In 11 calculations of the transverse cooling including multiple scatter- 
ing have also been performed. However, in order to present a fully self- 
consistent model, the authors have chosen to model the electron density by 
solving the Poisson equation for the harmonic potential. This yields an elec- 
tron density which is constant as a function of transverse position and may 
underestimate the net cooling effect as a cause of this. On the other hand, 
in 11 scattering 'on fluctuations of the planar potential', i.e. the nuclear 
contribution is not taken into account either. The analysis presented here 
supplements that of Baier and Katkov 11 in producing essentially the same 
conclusion by use of a different approach: Only under special circumstances 
may a penetrating particle experience a net cooling effect. 

5. Conclusions 

Qualitative agreement among two different theoretical approaches and an 
experiment on radiative angular cooling is shown. This gives considerable 
confidence in the predictions, especially of the semi-classical channeling ra- 
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diation reaction approach. Furthermore, for this approach multiple scatter- 
ing is included to obtain a more realistic estimate of the cooling properties. 
Rough estimates show that for existing positron beams a strong cooling 
effect may be achievable by means of axial channeling in a ~ 20 mm thick 
diamond crystal. There are however still open questions which deserve 
attention, among others that of dilution of the longitudinal phase-space 
due to the severe straggling in energy loss and the possibility of cooling 
simultaneously in both transverse directions by proper axial channeling or 
channeling in carbon nanotubcs 12 . 
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